Vegetation against dune mobility 
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Vegetation is the most common and most reliable stabilizer of loose soil or sand. This ancient technique is for 
the first time cast into a set of equations of motion describing the competition between aeolian sand transport and 
vegetation growth. Our set of equations is then applied to study quantitatively the transition between barchans 
and parabolic dunes driven by the dimensionless fixation index 8 which is the ratio between dune characteristic 
erosion rate and vegetation growth velocity. We find a fixation index 9 C below which the dunes are stabilized 
characterized by scaling laws. 

PACS numbers: 45.70.Qj,05.65.+b,92.10.Wa,92.40.Gc,92.60.Gn 
Keywords: 



Since pharaonic times mobile sand has been stabilized 
through plantations while conversely fields and woodland 
have been devastated by wind-driven erosion and coverage of 
sand. This ancient fight between vegetation growth and ae- 
olian surface mobility has evidently enormous impact on the 
economy of semi-arid regions, on coastal management and 
on global ecosystems. While empirical techniques have been 
systematically improved since dune fields in Aquitaine were 
fixed during the reign of Louis XIV and an entire specialty 
has developed in agronomy QJ 00 most approaches are just 
based on trial and error and there is still an astonishing lack 
of mathematical description and of quantitative predictability. 
It is the aim of this Letter to propose for the first time a set of 
differential equations of motion describing the aeolian trans- 
port on vegetated granular surfaces including the growth and 
destruction of the plants. Vegetation hinders sand mobility but 
also becomes its victim through erosion of the roots and cov- 
erage by sand. We will also show that these new equations 
can be used to calculate the morphological transition between 
barchans and parabolic dunes. 

Migrating crescent dunes occupy about half of the total 
desert and coastal area under uni-modal wind condition U 
and they exist either as isolated barchan dunes in places with 
sparse sand or forming barchanoid ridges where sand is abun- 
dant (Fig. 1, top left). These migrating dunes can be deac- 
tivated by the invasion of vegetation, as seen in Fig. 1 (top), 
provided that there is a certain amount of rainfall and weak 
human interference (HQ]. When the vegetation cover grows, 
barchans apparentely undergo a transformation into parabolic 
dunes with arms pointing upwind and partly colonized by 
plants (Fig. 1, top center-right), a metamorphosis that has been 
seen as the first step of dune inactivation |g,000]. 

The complex coupling between vegetation and aeolian sand 
transport involves two different time scales related, on one 
hand, to vegetation growth and erosion and deposition pro- 
cesses that change the surface, and, on the other hand, to sand 
transport and wind flow. A significant change in vegetation or 
in sand surface can happen within some hours or even days. 
In contrast, the time scale of wind flow changes and salta- 
tion process is of the order of seconds and therefore orders 
of magnitude faster. This separation of time scales leads to 
an enormous simplification because it decouples the different 
processes. Hence, we can use stationary solutions for the wind 



surface shear stress r and for the resulting sand flux q. 

The evolution of the sand surface and the vegetation that 
grows over it remains however coupled. Since plants locally 
slow down the wind, they can inhibit sand erosion as well 
as enhance sand accretion. This local slowdown of the wind 
shear stress represents the main dynamical effect of the veg- 
etation over the wind field and hence on the sand flux. The 
vegetation acts as a roughness element that absorbs part of the 
momentum transferred to the soil by the wind. As a result, 
the total surface shear stress r = |r| can be divided into two 
components, t v acting on the vegetation and r s on the sand 
grains. When plants are randomly distributed and the effec- 
tive shelter area for one plant is assumed to be proportional 
to its frontal area, the absorbed shear stress t v is proportional 
to the ve geta tion frontal area density A times the undisturbed 
shear r s jlQIl - Therefore, the fraction t s of total stress acting 
on sand grains is reduced to 

T S = r/(l + ml3X) , (1) 

where (3 is the ratio of plant to surface drag coefficients and 
A = p v /o-, where p v is the vegetation cover density and a the 
ratio of plant basal to frontal area, m is a model parameter 
flflfTlll . If the plant frontal area density A is zero there is no 
shear stress reduction. Otherwise, the reduction depends on 
the frontal area and on the plant drag coefficient encoded in (3. 

The time evolution of the sand surface height h is calculated 
using the conservation of mass, 

dh/dt = -V • q(r s ) , (2) 

defining the erosion rate dh/dt in terms of the sand flux q 
resulting from a wind surface shear stress r s that includes the 
vegetation feedback. 

Additionally, the sand dynamics affects the growth of 
plants, since non-cohesive sand is severely eroded by strong 
winds denuding the roots and increasing the evaporation from 
deep layers. Although desert plants can resist quite severe 
conditions, sand erosion often kills them I12l UJL fl4ll . We as- 
sume therefore that plants under sand erosion or deposition 
need extra time to adapt to the surface change and the vegeta- 
tion growth rate is delayed by the net erosion rate. Otherwise, 
they can grow up to a maximum height H v during a charac- 



2 




FIG. 1: Deactivation of migrating dunes under the influence of vegetation. On top, a dune field in White Sand, New Mexico, that shows 
barchanoid ridges on the left, where vegetation is absent, developing towards a mixture of active and inactive parabolic dunes on the right 
(wind blows from left to right). Dark green regions indicate abundance of vegetation. This suggests a transition between both types of dunes 
when the vegetation cover increases. This transition is illustrated with various dune types found in the White Sand dune field (pictures in 
the middle), reinforcing the idea of their common evolution from a crescent dune. Satellite images taken from GoogleEarth. Below, same 
transition obtained by the numerical solution of a model that accounts for the coupling between sand transport and vegetation, with fixation 
index 6 = 0.22. The vegetation cover is represented in green (grey). 



teristic growth time t v fl5ll 

dhy/dt = V v (1 - h v /H v ) - \dh/dt\ (3) 

where h v is the actual plant height, related with the vegeta- 
tion cover density introduced in Eq. \l]by p v = (h v /H v ) 2 . 
Besides, the vegetation growth velocity V v = H v /t v contains 
the information of climatic and local conditions that enhance 
or inhibit the growth process 1 1 2ll 1 3l 1 1 411 . 

This model for the interaction between plants and moving 
sand is closed by using an explicit formulation for the sand 
transport over a complex terrain 1 1 6ll 1 71 1 1 811 . Analytical cal- 
culations of the flow over a gentle hill yield an analytical ex- 
pression for the topographically induc ed p erturbation St of 
the surface shear stress f*o on a flat bed 11911 . therefore 

T(h) =f Q +T Q Sf(h), (4) 

where To = r*o|. Whereas the stationary sand flux q is given 
by the logistic equation 

V-q = q(l-q/q s )/l s , (5) 

where q = \q\, q s = \q s \, g s (r) is the saturated sand flux, 
i.e. the maximum sand flux carried by the wind with a sur- 
face shear r, and l s (r) is the characteristic length of saturation 
transients, called saturation length. 

Given an initial sand surface, a surface shear stress To on 
a flat bed and a sand influx, the time evolution of the com- 
plete system is determined by the consecutive integration of 
the equations for r(h) (Eq.gJ, t s (Eq.[TJ, q (Eq.0, h (Eq.|3 
and h v (Eq.[3}. 

Figures la-le show the effect of growing vegetation on 
the mobility of a barchan that evolves under a constant uni- 
directional wind (Fig. la). From Eq. [3] plants will grow 



wherever V v > \dh/dt\ even if under erosion conditions 
(dh/dt < 0) they can not germinate due to the mobile sand 
surface. Therefore, we introduce an heterogeneous seeding by 
allowing plants to germinate only where no erosion occurs, 
establishing a competition for their survival with the mobile 
sand. Such a scenario could for instance be the consequence 
of the cessation of human activity 0], an increase of the an- 
nual precipitation or a reduction of the wind strength in a dune 
field 171 12011 . all of them stimulating conditions for vegetation 
growth. 

Plants first invade locations with small enough erosion or 
deposition, like the horns, the crest and the surroundings of 
the dune. This vegetation reduces the strength of the wind 
and traps the sand which accumulates mainly where sand flux 
is highest, i.e. at the crest. Consequently, the sand cannot 
reach the lee side anymore and deposits on the crest (Fig. lb). 
This deposition kills the vegetation again. Thus the ero- 
sion/deposition rate at the crest is again lowered. This pro- 
vides a cyclic mechanism for the vegetation growth on the 
crest. 

As a result, the central part of the dune moves forward and 
two marginal ridges are left behind at the horns (Fig. lc) in a 
process that leads to the stretching of the windward side and 
the formation of a parabolic dune (Fig. Id). Finally, vegeta- 
tion overcomes sand erosion even in the central part and the 
migration velocity of the parabolic dune dramatically drops 
about 80% indicating its inactivation (Fig. le). This process 
agrees well with a recent conceptual model based on field ob- 
servations 0| . 

The transformation of barchans into parabolic dunes is de- 
termined by the initial barchan volume V, the undisturbed 
saturated sand flux Q = <f s (Yo)|, which encodes the wind 
strength in the flat bed, and the vegetation growth rate V v . 
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In fact, the stabilization process depends on the competition 
between sand transport and vegetation growth expressed by 
Eq.ljSJi- This competition can be quantified by a dimensionless 
control parameter 9, which we call fixation index, defined as 
the ratio between the characteristic erosion rate on the initial 
barchan dune, which scales as Q to the volume defined char- 
acteristic length V 1 ' 3 ratio (Eq.[2j, and the vegetation growth 
rate, 

9 = Q/iV^Vv) . (6) 

After performing simulations for different V, Q and V v , we 
find a common fixation index 9 C « 0.5 beyond which vege- 
tation fails to complete the inversion of the barchan dune. As 
long as 9 > 9 C barchan dunes remain mobile, otherwise they 
are deactivated by plants. 
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FIG. 2: Evolution of the normalized dune velocity (v/vo) for simu- 
lations with initial condition given by the parameters V (xlO 3 m 3 ), 
Q (m 2 /yr) and V v (m/yr): (1) 200, 286 and 78, (2) 200, 286 and 26, 
(3) 44, 529 and 52, and (4) 200, 529 and 26. They represent fixation 
indexes 6: 0.06 (1), 0.19 (2), 0.30 (3) and 0.35 (4), respectively. In- 
set (a): the normalized dune inactivation time t s /t m diverges when 
—* 9 C (solid line). Inset (b): the normalized final parabolic dune 
length L/V 1/l3 is proportional to t 3 /t m . Dots indicate simulation 
results. 

This morphologic transition at 9 C is typified by the evolu- 
tion of the normalized velocity v/vo of vegetated dunes for 
different initial conditions (some of them are shown in Fig. 2), 
where vq is the initial barchan velocity. They reveal a general 
pattern characterized for medium and high fixation indexes by 
a speedup preceeding a sharp decrease of the dune velocity at 
t = t s when inactivation occurs (the explicit condition for in- 
activation reads v/vq < 0.2. At t < t s the barchan evolves 
toward an active parabolic dune (Fig. 1, (a)-(d)) that becomes 
inactive for larger times (t > t s ) (Fig. 1, (e)). However, for 
fixation index larger than 9 C barchans remain mobile and thus 
t s diverges. 

Since the dune velocity v is proportional to the inverse of its 
height EUl . the sharp speedup in Fig. 2 is a consequence of the 
height decrease due to the redistribution of sand in the elon- 
gated dune. This redistribution process is enhanced at high 9 



when vegetation is progressively weaker against sand erosion. 
At small fixation indexes however, vegetation overcomes sand 
deposition and dune movement is determined by the vegeta- 
tion dying rate on the crest which results in a slower dune 
motion. In this regime the dune velocity oscillates due to the 
cyclic mechanism that controls the vegetation growth on the 
crest. The same mechanism is responsable for the short reac- 
tivation of the inactive parabolic dune represented by the oc- 
casional second peak in the dune velocity for t > t s (Fig. 2). 

From the analysis of the inactivation time t s as function of 
V, Q and V v , we find that t a scales with a characteristic dune 
migration time t m = V 2 / 3 /Q and obeys a power law in 9 
(Fig. 2, inset(a)) 

t,»0.17i m /(e c -9), (7) 

which discloses a transition between an inactive parabolic 
dune and a barchan dune remaining active. 




FIG. 3: Examples of parabolic dune shapes at fixation indexes 9; (a) 
0.16, (b) 0.22 and (c) 0.27. (d) is an example of an intermediate state 
at 8 = 0.38 where the dune still is barchan-like, in sharp contrast 
to the rest. Each figure is compared with real examples from White 
Sand, New Mexico. Note the similarities in the contour lines market 
in red. Satellite images taken from GoogleEarth. 

The normalized final length of the parabolic dune L/V 1 ^ 3 
also diverges near the transition point. The parabolic 
length depends linearly on the normalized inactivation time, 
L/V 1/3 w at s /t m - r, (Fig. 2, inset (b)) where a w 24.5 
and r ~ 6.2, and thus scales as L/V 1/3 ~ (0 C - 9)" 1 . Note 
that in the limit 9 = 0, the inactivation time is still propor- 
tional to t m (t s ?s 0.34 t m ) and the dune moves a distance 
L k 2 V 1 / 3 . This is consequence of the heterogeneous seed- 
ing assumed before. Since plants can not grow in the dune 
windward side where erosion occurs, all this sand is stabilized 
only after crossing the crest. Thus, the dune moves a distance 
equivalent to its windward side length ~ 2 V 1 ! 3 . 

Figure 3 shows some examples of the L dependence on 9. 
For small 9, barchan dunes are quickly deactivated (EqQ and 
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short parabolic dunes emerge (Fig. 3a). At large 9 the amount 
of sand trapped by plants in the dune arms, which height 
scales as V 2 / 3 /L ~ V x l 3 (9 C — 9), is reduced. This leads 
to longer parabolic dunes whose 'noses' can experience suc- 
cessive splits before being finally deactivated (Fig. 3b and 3c). 
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FIG. 4: The fraction of available sand surface 1 — (p v ), where (p v ) 
is the mean vegetation cover density, shows an exponential decay 
during the evolution of a vegetated barchan dune for small 8 when 
time is normalized by t m . Inset: proportionality between the mean 
sand flux (q), normalized by the saturated upwind value Q, and the 
fraction of active sand surface. 

The dependence of t s on the properties of the initial barchan 
through 9 and t m leads us to characterize the mobility or activ- 
ity of a vegetated dune only based on general field conditions. 
In particular, the condition 9 < 9 C for a migrating barchan to 
become inactive explicitely reads 

Q <Q C = 9 C V V V 1 / 3 , (8) 

which represents a size dependent upper limit for Q below 
which inactivation takes place, and thus stresses the lead- 



ing role of the wind strength in the inactivation process. It 
also defines a minimum vegetation growth velocity V vc = 
9~ 1 Q/V 1 / 3 , also a function on dune size, and dune volume 
V c = (9~ 1 Q /V v ) 3 above which barchan dunes are stabilized. 
Furthermore, the dune inactivation time as explicit function of 
V achieves a minimum for a dune volume Vq ~ (3/2) 3 V c . 
Hence, dunes with volumes below V c will remain mobile, 
while dunes with volumes Vq just above the minimum one 
will be the first to be stabilized. Large dunes however, will 
remain mobile for a time that scales with t m oc V 2 ! 3 . 

The dynamics of dune stabilization is also characterized 
by the evolution of the mean vegetation cover density (p v ) 
(Fig. 4) and, in a minor way, by the mean sand flux (q) over 
the dune, which we find to be proportional to the fraction of 
available sand surface 1 — (p v ) (Fig. 4, inset). For small 9 
the evolution follows the same trend, e.g. a rising vegetation 
cover and the resulting reduction of the mean sand flux I22I1 . 
with t m as characteristic time. For larger 9, the vegetation 
cover reaches a plateau that extends until the dune stabilizes 
at t ~ t s , after which vegetation cover rises again due to the 
inactivation of the dune. This plateau extends indefinitely at 
9 = 9 C , when t s diverges, as a sign of actual sand mobility. 

We have proposed a set of differential equations of mo- 
tion describing the aeolian transport on vegetated granular sur- 
faces including the growth and destruction of plants. Through 
them we calculate the morphological transition between active 
barchans and inactive parabolic dunes characterized by the di- 
vergence of the dune inactivation time t s and its scaling with 
the barchan migration time t m . We find that the fixation index 
9, a dimensionless combination of barchan dune size, wind 
strength and vegetation growth velocity, determines from the 
very beginning the final outcome of the competition between 
vegetation growth and aeolian surface mobility. These predic- 
tions well be directly tested in the field and have implications 
on the economy of semi-arid regions on coastal management 
and on global ecosystems. 
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